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1. INTRODUCTION 

The idea that planetary atmospheres can erode as a result of impact, and thus lose 
mass along with solid and molten high velocity ejecta during accretional infall of 
planetesimals follows from such early thoughtful works as that of Arrhenius et al [1974], 
Benlow and Meadows [1977], Ringwood [1979], and Cameron [1983]. Ahrens et al 
[1989] describe how planetary impact accretion (and impact erosion) concepts lead 
naturally, from the idea that atmospheres form and erode, during planetary growth. 

The theory of planetary system formation from a disc of gas, and later, gas and 
dust, corotating around the proto-sun, is described by Safronov [1969], Wetherill [1980], 
and Kaula [1979] (Fig. 1 ). In this model, planets grow as a result of mutual attraction 
and collision of planetesimals within the early solar nebula. An important step in the wide 
acceptance of this model was the theoretical work of Goldreich & Ward [1973] who 
quantitatively showed that condensed material from a solar nebula would grow into the km- 
planetesimals hypothesized by Safronov. Recently, very strong support for the 
Safronov- Wetherill-Kaula scenario of planetary growth via mutual collision has come from 
discovery using infrared, optical, and radio imagery of disc-shaped circumstellar gas and 
dust clouds around T-Tauri stars (e.g. HL Tauri and R Monocerotis) [Beckwith et al, 1986; 
Sargent & Beckwith, 1987] and main sequence stars (e.g. jJ-Pictoris) [Smith & Terrile, 
1984] (Fig. 2). Although the gaseous rotating discs around T-Tauri stars have been 
imaged via microwave radio interferometry and infrared techniques, it has not been 
demonstrated that these discs contain accreting planetesimals. Recent observations of the 
13co emission spectra and the blackbody thermal emission from the circumstellar gas disc 
around HL Tauri indicate that this disc extends out to a radius of 2000 AU, but is less than 
380 AU thick. Moreover, the spectra are consistent with the gas and dust moving in bound 
orbits around this star [Sargent & Beckwith, 1992]. 

I review current ideas about the nature of the planetesimals —composition, size 
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distribution, and the planetary encounter velocity. Previous papers on accretion and 
erosion of planetary atmospheres as a result of multiple impacts are reviewed. Finally, the 
effects of blowing off a substantial fraction of the atmosphere from a terrestrial planet due 
to a single giant body impact are discussed. 

2. PLANET FORMING MATERIALS 

The planets and minor objects in the solar system appear to have accreted from the 
following three components: 

1) Planetesimals similar to meteorites. 

The constitution of the terrestrial planets suggest that they accreted largely from 
planetesimals with a range of composition including primitive objects, such as. Cl 
chondrites, as well as objects similar to differentiated metal and silicate meteorites. The 
linear relation in Fig. 3 indicates how similar the major element composition of the 
terrestrial planets are to the sun, while Fig. 4 illustrates that the noble gas abundance 
patterns of terrestrial planets are similar to those of primitive meteorites such as Cl 
chondrites. Meteorites are often taken to be typical of the planetesimals existing within the 
inner zone of the solar nebula from which the terrestrial planets accreted. The vestiges of 
the planetesimals of the inner solar system are believed to be the asteroids. We presumably 
sample these objects via meteorites that fall on the Earth. Some of the planetesimals which 
formed the terrestrial planets were also probably similar in composition to the present 
comets. Largely because of the gravitational perturbation due to the large masses of Jupiter 
and Saturn, comets can achieve highly eccentric orbits which, near perihelion, result in 
collision with the terrestrial planets. Shoemaker et al [1990] estimated that some 30% of 
the recent impactor flux on the earth may be cometary. How much of this flux has 
provided the volatile budgets of the terrestrial planets is presently unclear (e.g. Grinspoon 
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& Lewis [1988] and Donahue & Hodges [1992]) because of the similarity of the cometary 
non-volatile major element inventory with that of the Cl chondrites [Jessberger et al, 

1989 ]. 

2) Cometary Planetesimals 

In contrast to the silicate and iron-rich planetesimals, the planetesimals that made up 
of the cores of the giant outer planets — Jupiter and Saturn and to a lesser degree, Neptune 
and Uranus, and possibly, Pluto, were ice-rich. Their remnants are now associated with 
comets. The generally smaller size (< 10 km diameter) of the comets (relative to asteroids) 
suggest that in the outer solar system the density of matter in the solar nebula was never 
great enough for the ice-rich planetesimals to experience substantial mutual gravitational 
attraction, resulting in impact accretion. 

3) Solar Nebulae Gas 

In contrast to terrestrial planets, Jupiter and Saturn, and to a lesser degree, Uranus 
and Neptune appear to have, after building up their initial planetary core from ice, silicate, 
and carbon-rich objects, gravitationally captured large quantities of solar nebular gases. 

The large planets thus retain a large solar-like reservoir of H2 and He. Pollack & 
Bodenheimer [1989] suggested that the ratio of carbon to hydrogen in the atmospheres 
which dramatically increases in the planets, in the order: Jupiter, Saturn, Uranus, and 
Neptune, directly reflects a decreasing budget of accreted solar nebular hydrogen. 

3. TERRESTRIAL PLANETARY VOLATILES 

On the basis of the general similarity of the noble gas abundances of the terrestrial 
planet atmospheres to the noble gas component which was processed within planetesimals 
before their accretion, I infer that terrestrial planet atmospheres originated from 
planetesimals similar to primitive meteorites. Cometary-like planetesimals may have also 
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contributed to the terrestrial planet volatile inventory. Notably, the thermal and 
gravitational evolution of the planets are distinctly different, however, their relative noble 
gas inventories are similar. The solid phases containing noble gases presumably are still 
present in the asteroids and demonstratively occur in meteorites. Small meteorites impact at 
sufficiently slow terminal velocity that the noble gases are not released upon impact They 
still contain the complement of noble gases which were present 4.6 Gyr ago prior to the 
planet-forming epoch (Fig. 2). Moreover, the two key mass-selective, gas-loss 
mechanisms (e.g. Hunten et al [1989]) for atmospheric escape from planets, Jeans loss and 
hydrodynamic escape, predict abundance patterns, starting from a solar noble gas pattern, 
which are quite different from those in the atmospheres of Mars, Earth and Venus. 

4. SIZE DISTRIBUTION OF PLANETESIMALS 

In Safronov's [1969] theory of planetary accretion, a simple power law distribution 
of the differential number of objects, dN, which occurred in a mass range, dm is assumed 

dN sm-9 dm (1) 

where 1 < q < 2. 

For veiy energetic collisions, more mass becomes concentrated in the fines, which 
implies that q > 2. As pointed out by Greenberg [1989] when q = 2, each size range 
contains exactly the same total mass. Observations of the present asteroids, which 
presumably have collided and reaccreted over their 4.6 Gyr lifetime show a mass 
distribution close to q = 2. However, statistical studies of accretion models indicate that 
values closer to q = 1.5 occurred in the past which implies that most of the mass of the 
impactor population occurred in the largest plane tesimals. Weak constraints on the value of 
q for impactors in the solar system come from the distribution of crater sizes on ancient 
cratered terranes, such as the highlands of the Moon. If a planet with a current mass, M, is 
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accreting a veneer of material, of mass, 8 M, the mass of the largest impactor, mi, is 
obtained by integrating: 


8M = [ m, mN(m)dm (2) 

and assuming that Eq. 1 is valid up to the size of the largest bodies, for q < 2 [Greenberg, 

1989; Wetherill, 1985; Wetherill & Stewart, 1989] 

mj - 8M (4 — 2q) / (3 — q) (2a) 


Thus, for example, if the accretion of 10% of the earth, 8M = 6x10^6 g and q = 1.5, then 
mi = 4X10 26 g, or mi has a radius of 2.6 x 10 8 cm, which is similar to Mar’s radius. Thus 
we need to consider the effect of impact on planets of very substantial projectiles having 
radii, perhaps 10 2 times the scale height of the planetary atmosphere. 

5. PLANETESIMAL IMPACT VELOCITIES 

As shown in Fig. 1, the major element of motion of the planetesimals in ideal 
circular Keplerian circular orbits, around the sun, is velocity, v$, which, in the solar 
reference frame is given by 

v s = ^GM e / r (3) 

where G is the gravitational constant (6.67 x 10* 8 dyn-cm 2 /g 2 or 6.7 x 10' 11 nm 2 /kg 2 ), M© 
is the solar mass and r is the distance of the planetesimal from the sun's center. Even for 
perfectly circular planar orbits, two objects at slightly different solar distances, initially will 
differ in velocity. Hence, the object closest to the sun, mj, will move in at a slightly higher 
speed, passing an object, m2, which is further from the sun. Gravitational interaction 
occurs in this two-body encounter, such that the object, mi, will experience a radial 


velocity increase given by 
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8v = m 2 vj/ (mi + m 2 )E (4) 

Here, E, is an encounter parameter which depends on the geometry and the relative masses 
of the objects. Continual gravitational interaction of adjacent objects gives rise to 
increasing orbital eccentricities, as well as inclination of the orbits. These in turn give rise 
to collisions which tend to damp out the velocity and orbital perturbations such as described 
by Eqn. 4. On average, for planetesimals relative to one another in orbit in a corotating 
disc of particles, their encounter velocities will be v«>. Safronov pointed out that the largest 
particles with an escape velocity, v e , will have encounter velocities: 

Voc = Ve/V 2 0 (5) 

Here the local escape velocity is 

v c = ^2Gm p /R p = J2 R p g (6) 

The result depends critically on q < 2. Otherwise, v« will be drastically reduced by 
multiple encounters with the smallest particles. 

The impact velocity of a planetesimal is therefore: 

Vi = V(v e 2 + v« 2 ) (7) 

where m p is the mass of the largest planetesimal in a region, R p is the largest planetesimal 
radius, and g is its gravitational acceleration. Here, 0, the Safronov parameter, is usually 
taken to be about 4 or 5. Eq. 7 implies that as planets grow by accretion, the velocity is 
always somewhat greater than the planetary escape velocity. Moreover, Eq. 2a indicates 
that as planets grow, so does the mass of the planetesimals which impact their surface. 
These considerations have led to efforts to understand the essential physics of large body 
impacts on the terrestrial planets by Benz et al [1989; 1986; 1987; 1988] and Kipp & 
Melosh [1986]. 
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6. COACCREnON OF PLANETARY ATMOSPHERES 

Although Lange & Ahrens [1982b] suggested that the impact-induced dehydration 
of water-bearing minerals in planetesimals such as serpentine Mg3 S^Qs (OH)4 would 
produce a largely water-rich atmosphere on the growing planets, it was Abe & Matsui 
[1985] who first suggested the possibility that water in this atmosphere, and possibly the 
dust produced by planetesimal impact could drastically alter the thermal regime on the 
surface of growing planets (Fig. 5). They assumed that serpentine in planetesimals 
brought the Earth and the other terrestrial planets their water inventory during their 
accretion. Once impact pressures of infalling planetesimals exceeded P ~ 23 GPa (for a 
porous regolith), the supply of water to the atmosphere was assumed to begin. Abe & 
Matsui [1985] estimated the peak shock pressure by 

P = Po [Q> + (K' + 1) Vj/8 ] Vj/2 (8) 

where po and Co are impacting planetesimal initial density and zero-presssure bulk sound 
velocity, respectively, and K' is the partial derivative of the bulk modulus with respect to 
pressure. In these atmospheric accretion models the shock pressure is assumed to act on 
the entire impacting planetesimal and H2O is released, as well as lesser amounts of CO2, 
NH3, SO2, and other volatiles present in primitive meteorites (e.g. Cl carbonaceous 
chondrites). Previous calculations have also included contributions to the atmosphere 
induced by shock-loading of volatiles already present in the material of the planetary 
surface layer (e.g., Lange & Ahrens [1982a]). 

Using the experimental value of 23 GPa as the shock pressure required to induce 
complete water loss for serpentine Lange & Ahrens [1982b] and Abe & Matsui [1986] 
concluded that once the radius of the Earth reached between 0.2 to 0.4 of the present value 
(Re), thermal blanketing of the Earth caused by a dense water atmosphere occurred. 

Matsui & Abe [1986a] point out that thermal blanketing is a more severe condition than the 
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greenhouse effect In a greenhouse effect sunlight penetrates through the atmosphere (in 
die visible) but the thermal energy to be reradiated by the planetary surface in the infrared is 
trapped because of the infrared opacity of the planetary atmosphere resulting from abundant 
CO 2 and H 2 O. Thermal blanketing is more severe, because solar radiation incident on the 
top of the atmosphere is completely scattered by the fine aerosols and impact ejecta and all 
the thermal energy is absorbed by the greenhouse gases, H 2 O and CO 2 . Impact cratering 
calculations [Ahrens et al, 1989; O'Keefe & Ahrens, 1977a] demonstrate that for large 
impactors, 60 to 90% of the energy of the impact is delivered as internal energy of the near 
surface material. The major effect of the proto-atmosphere then is to provide an insulating 
blanket to the flux of heat due to impacts on the surface of the growing planet These 
processes are described by the equation 

Gm n 

(1 + 1/26) -mdt = 4nR 2 (F atm - F,)dt 

R p 

+ Cpriip (T s — Tp) dt + C p m s T s dt (9) 

where the left-hand side is the rate of kinetic energy supplied to the surface provided by the 
impacting planetesimals. The first term on the right is the balance of heat flux of the 
atmosphere where Fatm is the energy flux escaping from the surface to interplanetary space 
and Fi is the energy flux from the interior to the surface layer. The second and third terms 
are the heat sinks to the planet as a result of heating a larger planetary material of mass m s 
from planetesimal temperature T p to the higher surface temperature T s . For a graybody 
radiative equilibrium atmosphere 

Fatm = 2«»T« 4 - S 0 /4)/[(3kp 0 /2g) + 2 ] (10) 

where a is the Stefan-Boltzmann constant. So is the solar flux, k is the absorption 
coefficient in the atmosphere, and p 0 , is the surface pressure. Thermal blanketing as 
described by Eqs. (9) and (10), with the reasonable values of the constants chosen by Abe 
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& Matsui [1985; 1986], quickly leads to temperatures above the solidus of crustal (basaltic) 
rocks (Fig. 6a). Moreover, as more planetesimals impact the planet, the additional water 
provided begins to dissolve in what is the start of a magma ocean. Abe & Matsui showed 
that the surface temperature should be buffered by the solidus of hydrous basalt, - 1500 K, 
and the mass of H 2 O in the atmosphere is nearly constant at 10 24 g (essentially the present 
surface water budget) (Fig. 6b). They showed that the mass of H 2 O in the atmosphere 
remains nearly constant (due to a negative feedback effect). This effect can be 
demonstrated for a small surface temperature increase. This causes an increase in the 
fraction of molten basalt That in turn induces additional water to dissolve in the molten 
silicate. Loss of water from the atmosphere then decreases the effectiveness of atmospheric 
blanketing and the result is that the small increase in temperature is nullified by the system's 
negative feedback. 

The effectiveness of thermal blanketing of impact energy by the massive 
protoatmosphere, as well as the feedback effect of water solubility in molten silicates has 
been independently verified in a study by Zahnle et al [1988]. 

The termination of the coaccretion of an atmosphere and planet, which was modeled 
to occur on Earth, Venus and Mars, can occur via three different mechanisms. 

1) Abe & Matsui [1988] suggested that as the accretion rate decreased, the impact 
energy flux at the base of the atmosphere decreased and gradually solar heating dominated 
over impact heating. The surface temperature then declined below the melting point of 
hydrous basalt With decreasing temperature, the water condensed and formed terrestrial 
oceans. Oceans may have formed on Venus [Matsui & Abe, 1986b], however, the larger 
solar ultraviolet flux gives rise to a enhanced photodisassociation to hydrogen and oxygen 
in the upper atmosphere and subsequent Jean's escape of a large fraction of the planet's 
hydrogen inventory gives rise to the presently observed enhancement of the D/H ratio of 
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Venus relative to the Earth of -10 2 . 

2) Atmospheric loss occurs via multiple impact erosion. This is discussed in the 
next section. 

3) Sudden partial or complete atmospheric loss occurs as a result of a large body 
impact This is discussed in Section 7. 

7. ATMOSPHERIC EROSION BY IMPACT CRATERING 

In addition to bringing volatiles to accreting atmospheres, the infall of planetesimals 
can erode planets and their oceans and atmospheres. For atmosphere-free solid and molten 
silicate planets, because the mechanical impedance of the impacting planetesimals is, in 
general, similar to that of a planet the amount of ejecta which can escape from a planet with 
a given surface escape velocity depends only on impact velocity (and hence energy per unit 
mass) [O’Keefe & Ahrens, 1977b], whereas the net gain or loss of a planetary atmosphere 
depends on total impact energy. O’Keefe & Ahrens [1982] calculated the energy 
partitioning into an atmosphere overlying a planet They found that upon impact of a planet 
with projectiles of radii less than the atmospheric scale height where the ejecta was 
primarily solid or molten, the amount of energy imparted to the atmosphere by direct 
passage through the atmosphere was only a few percent Moreover, very little of the 
atmosphere achieved upward velocities in excess of the escape velocity. Walker [1986] 
showed that a very small portion of the atmosphere, shocked by the meteoroid, achieved 
sufficient enthalpy density to expand to greater than escape velocity. Using a numerical 
explosion model, Jones & Kodis [1982] showed that for the Earth, atmospheric explosion 
energies > 5 x lO 2 ^ ergs induced significant atmospheric blow-off. 

Subsequently, Ahrens & O'Keefe [1987] and Ahrens et al [1989] employed a 
model in which they assumed all the energy of the impactor is delivered to the planetary 



13 


surface and applied a theory for the shock acceleration of the atmosphere by an explosion 
developed by Zel'dovich & Raizer [1966 ChapL 12] and Bach et al [1975] (Fig. 7). In this 
model, the time for atmospheric escape is related to atmospheric density near the Earth’s 
surface, poo, explosion energy, E, and atmospheric scale height, H, by 

t * Ci (poo H 5 /E) 1/2 (11) 

where the constant, Cj is approximately equal to 25. Moreover, the initial atmospheric 
shock velocity for atmospheric escape is 

D = aH/t (12) 

where a = 6. Assuming the Earth impactor has a velocity, v e = 11 km/sec and a strong 
shock condition exists such that D = 1 1 km/sec, Eq. 12 yields t = 4.4 sec and the minimum 
impact energy calculated from Eq. 1 1 is 1.9 x lO 77 ergs. A projectile carrying this energy, 
if composed of silicate will have a radius of - 0.5 km. The energy, 1.9 x 1(£ 7 ergs, is 
somewhat greater (by a factor of 20) than from the numerical calculations of Jones & Kodis 
[1982]. 

For more energetic impacts, (larger impactors) Melosh & Vickery [1989] and 
Vickery & Melosh [1990] developed a simple atmospheric cratering model applicable, for 
example, on the Earth in the 4.5 x 1(£ 7 to 9.9 x lO 3 ^ erg energy range (Fig. 8). For 20 
km/sec, 2.7 g/cm 3 impactors, these energies correspond to impactor radii of 0.6 to 7.6 km. 
They pointed out that for high-velocity planetary impacts, which penetrate the atmosphere, 
the projectile and an equivalent mass of target become vaporized. The resulting gas plume 
then expands at a speed greater than the planetary escape velocity and carries with it the 
overlying planetary atmosphere. A conservative model of the plume expansion gives the 
mean (mass-averaged) velocity of expansion as 

1/2 


v exp “ [^( e hyap)] 


(13) 
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where e is the kinetic energy per unit mass of the impactor and hyq> is the total enthalpy per 
unit mass, starting at ambient temperature required to vaporize the projectile or target 
material. For silicate and ice, hyap is 13 x 10*3 and 3 x 10 1 * ergs/g, respectively. The 
energy per unit mass imparted by the impact-induced shock-wave for like materials is e - v 
2/8, where we assume the shock particle velocity, u is v/2. It is easy to show that the 
minimum impact velocity required for the vapor plume to exceed the escape velocity is 

v m = A /8( V 2/2 + h vap ) (14) 

Atmospheric erosion occurs if the impact-induced momentum of this shock-induced 
gas when combined with the mass of the overlying atmosphere has sufficient velocity to 
escape the planet By assuming a velocity profile in the total expanding gas cloud proposed 
by Zel’dovich & Raizer [1966 p. 104] of the form 

p(r) = A( l-r2/R2)ot/R3 (15) 

where R is the radius of the front of the gas cloud, r is the radius to a point within the gas 
cloud, and A and a are determined by assuming conservation of mass and energy above 
the impact site. Vickery & Melosh [1990] used a value of y= 9/7 to infer a value of a = 

1 1 and A = 15.4M t , where y is the gas cloud’s polytropic exponent, and M t is the total 
mass of the vapor cloud. When the projectile velocity exceeds the minimum impact 
velocity for atmospheric escape (Eq. 14), escape occurs in atmosphere directly above the 
impact site. As the impact energy is increased, a cone with increasing angle 8 is ejected 
(Fig. 9). As impacts become more energetic, the maximum energy in the Vickery-Melosh 
model corresponds to ejection of an air mass above the tangent plane (Fig. 9) of 3.3 x 10 19 
g (for the case of the Earth) or 6.3 x 1(H of the total atmospheric budget Vickery and 
Melosh showed that when the mass of the projectile exceeds the mass of the atmosphere 
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above the horizontal tangent, me, all of the atmosphere above a tangent plane to Earth is 
ejected. Thus when 

m ^ me * H m^2Rp 

Eq. 16 thus indicates that for smaller m a , the mass of the atmosphere, gives nse to a 
smaller mass, me, which can erode an atmosphere. Table 1 gives me for the terrestrial 
planets. Since smaller projectiles are more numerous and thinner atmospheres erode 
rapidly, Eq. 16 indicates that once an atmosphere starts eroding, erosion is accelerated until 
the planet is stripped. Fig. 10 shows the maximum atmospheric mass that can be expelled 
by a spectrum of impactors for three different planets for a given accretion (veneer) mass 
being added to the planet Melosh & Vickery [1989] demonstrated how Mars’ atmosphere 
could have been eroded from an initial surface pressure of 0.7 bars to the present 7x 10* 

3bars in the first 1.4 Gyr of solar system history. 

In conclusion, I note that in the case of the Earth, impact erosion is important only 
for impacts more energetic than -id* ergs. Thus the Abe-Matsui scenario, which only 
deals with atmospheric accretion, is relevant if the projectiles impacting the Earth’s surface 
axe much smaller in radius than the atmospheric scale height For projectiles with radii m 
the km range, the research summarized in this section indicates that impact erosion needs to 
be taken into account in accretion models of the terrestrial planets. 

8. ATMOSPHERIC BLOW-OFF BY GIANT IMPACTS 

As discussed in Sect 3, as the planets accreted according to the Safronov- 
Wetherill-Kaula scenario, the planetesimal impactors also grew in size and it appears likely 
that some planetesimals grew to radii in the 2,000 to 3,000 km range which are comparable 
to those of the smaller planets (e.g. Mercury and Mars) . These considerations have 
motivated the numerical modeling of large body impacts by Benz et al and Kipp and 
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Melosh cited in Section 4. 

Previously, the effects of large body impact on the Earth’s atmosphere has only 
been briefly described by Ahrens [1990]. To calculate the energy, and hence, approximate 
plane tesimal size, such that upon impact the entire planetary atmosphere is blown off, 1 
employ a different approach than previous efforts and consider a shock wave that is entirely 
propagated within a terrestrial planet as sketched in Fig. 11. 

The key calculation is to relate the particle velocity of the solid planet-atmosphere 
interface, Uf S , antipodal to a major impact, to the atmospheric free-surface velocity, v e , 
which will result from being shocked first by the moving solid planet and then 
isentropically being released into space at speeds greater than the escape velocity. This is a 
conservative calculation since I use the density and pressure of the atmosphere at its base. 
Moreover, the atmosphere covering the planet closer to the impact than the antipode is 
expected to achieve high velocity because it is shocked by the decaying air wave, and also 
shocked to higher pressures by the solid planet Note that, in general, as a shock wave is 
propagated upward in an exponential atmosphere, because the density encountered by the 
traveling shock is decreasing, the shock velocity and particle velocity increases with altitude 
as discussed in Sect 6 and by Zel'dovich & Raizer [1966 Vol. 2, Chapt XII, section 
5.25]. Thus, one can safely neglect shock attenuation in the atmosphere, and assume the 
particle velocity at the solid planet-atmosphere interface (the independent variable) and 
calculate the shock pressure induced in the gas by the outward surface of the Earth. The 
solid Earth therefore acts like a piston with velocity, Uf S , pushing on the atmosphere. 

The pressure behind the shock wave, pi, for different outward rock velocities, Uf S , 
can be determined by solving for pi in the following set of equations (e.g. eqs. 1.78+1.79, 


Z+R): 
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IUf s l = Uo-Ul ( 17 ) 

where 

% = {V 0 [(Y-l) Po + (Y + D Pll/2) 1/2 (18) 

is the particle velocity of the unshocked gas with respect to the shock wave in the 

atmosphere and 



v 0 [(Y + l) Po + (Y-l) Pi' 

2 ' 

U 1 = * 

2 [(Y-l) Po + (Y + l) Pi 



(19) 


where ui is the particle velocity of the shocked atmosphere relative to the shock front 
Surface values of p 0 and Vo=l/poo for the planetary atmospheres used to calculate pi are 
given in Table 2. Upon isentropic release from pi to zero pressure, an additional large 
increase in particle velocity occurs which is given by the Riemann integral. 

u, =f p '(-<JV/dp) 1 ' 2 dp (20) 

Jo 

Upon substituting for an ideal polytropic gas, with a polytropic exponent, y, Eq. (20) 
yields 

Ur = (Po 1 * Vo/y)i/2 pi i/2- 1/(2 y)/[ 1/2 - l/(2y)] (21) 

The shock-induced outward atmospheric velocity, %, is given by: 

Ue = hifcl + u r (22) 

We assume that when 

lue£v c l (23) 

atmospheric blow-off occurs. When Uc = v e , the corresponding value of u fe is denoted by 
Ufse- The outward rock velocity versus outward atmospheric velocity for the Earth, Venus, 
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and Mars, is shown in Fig. 12 for values of the polytropic exponent in the range from 1. 1 
to 1.3. This range encompasses the effective likely range of y which is expected to 
decrease from 1.3 to 1.1 with increasing gas ionization. Eq. 23 is satisfied for y - 1.1 to 
1.3 for outward rock velocities of 1.60 to 2.25 km/sec, 1.00 to 2.45 km/sec and 0.27 to 
1.2 km/sec for the Earth, Venus, and Mars. What impact energies will produce these 
outward rock velocities for the terrestrial planets? 

Fortunately, the strength of the shock-wave induced upon propagation completely 
through planets with varying iron core sizes overlain by silicate mantles has been recently 
calculated for objects which have core to planetary radius ratios of 0.333 and 0.466 [Watts 
et al, 1991]. In these calculations, the energy of the surface source was Ew = 3.1 xlO 34 
ergs. Notably, the calculation of Watts et al [1991], when scaled as discussed below, 
agrees closely with those of Hughes et al [1977] upon which earlier estimates of the energy 
required to blow-off the Earth's atmosphere were based [Ahrens, 1990]. 

I scaled the results, which are given as peak compressional wave stresses 
experienced by material directly beneath the antipode of the impact point for a core to 
planetary radius ratio of 0.333 and a planet radius, R w = 1500 km, to that for Mars, which 
has a core to radius ratio of -0.34 and a planetary radius of Rp = 3396 km (Table 2). 
Similarly, I used the Watts et al. result for a core to planetary radius ratio of 0.47 to provide 
estimates of the peak stress beneath the antipode for an impact on the Earth and Venus for 
which the actual core to planetary radius ratios are 0.44 and 0.47, respectively. Since 
planetary gravity was not included in the calculations, we employ cube scaling [Melosh, 
1989 p. 1 12] to adjust the results for planetary size. 

The energy of the equivalent surface source, Ep, assumed for an impact on the 


actual planet of interest is: 


E p = E w (Rp/R w )3 


(24) 
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To relate the peak pressure experienced by the cell beneath the antipodes for the 0.33 and 
0.47 core to planetary radius ratios, a peak shock pressure of Pj = 2.03 and Pi = 2.10 GPa 
were used. To convert these values to shock particle velocity, ui, I assumed a surface 
density of p 0 = 2.72 g/cm 3 and a shock velocity U s =5 km/sec in the momentum equation: 

ui = Pi / (Po U 8 ) (25) 

and then made the common approximation that the outward rock (free-surface) velocity 
Ufcw corresponding to the Watts et aL calculation is 

Ufsw = 2 ui (26) 

The energy, Efs required of an impactor to obtain the upper and lower bounds of ufs 
necessary to launch the atmosphere to escape velocity in the calculation of Fig. 12, can then 
be calculated from 

E&= Ep (ufc/ufcw) 2 ( 27 ) 

Thus, for complete atmospheric blow-off, values for Efs of 6.22 x 10 37 to 1.84 x 10 38 ergs 
are necessary for the Earth. This compares to -10 37 ergs previously calculated by Ahrens 
[1990]. For Venus, Ef S varies from 2.1 x 10 36 to 1.3 x 10 38 ergs, whereas for Mars 
values of 2.9 x 10 35 to 5.7 x 10 36 ergs are needed for complete blow-off (Table 1). 

Finally, it is useful to estimate the energy and mass fraction of the planetary 
atmosphere blown-off in going from the Melosh- Vickery model of tangential blow-off 
(Table 1) to the condition of complete blow-off. For energies less than those required to 
eject the entire atmosphere, I assume a simple power law for the decrease of particle 
velocity with radius, r, from the impact point (Fig. 1 lb) and assume that as the decaying 
stress wave interacts with the free-surface of a spherical planet, atmospheric blow-off 
occurs if the particle velocity is greater or equal to 0.5 uge, which is the value of Uf s , 
required to blow off the atmosphere for the antipodal case. Thus, the amount of 
atmosphere blown-off is related to the area of a sphere subtended by an arc of radius, r, at 
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the point where the particle velocity has decayed to 0.5 u&e (Fig. 1 lb). 

To determine an empirical relation for the attenuation of particle velocity, I assume 
the form [Melosh, 1989 p. 62]. 

u=Uo/r® (28) 

I first calculate the radius of a hemisphere, rj, enclosing a unit mass 

n =[3/(2* p 0 )] 1/3 (29) 

Using the values of Ef s min and Ef s ^ calculated from eq. 27 to designate the minimum 
and maximum energies obtained from Fig. 1 2 for y = 1.1 and y= 1.3, respectively, the 
shock particle velocity associated with each energy is given by 

U1 min = (Efs min)^ (30a) 

UJ max = (Efs max ) 1/2 (30b) 

Denoting luf S l calculated from Eq. 17 using y = 1.1 and y= 1.3 as uf s m m and Uf sma x I can 
obtain from Eqs. 26 and 28, expressions for the particle velocity decay parameters nj and 
n 2 from 

nj = log (2 uj min / Ufs min) / log (2 Rp/rj) (31a) 

n2 = log (2 uj max s max) / log (2 Rpfrj) (31b) 

For stress wave particle velocity decaying to a value of Uf sw /2 at a radius from the impact 
point r c , the mass of atmosphere blown off is 

me * Wg) A (32) 

where the term in parenthesis is the atmospheric mass per unit area and A is the area of a 
sphere subtended by an arc of length, r, (Fig. 1 lb). 

From geometrical arguments it can be shown that the area of the planet subtended 
by an arc of length, r, is 

A = 7t r 2 (33) 

The energies (minimum and maximum) associated with the radius, r c , for y = 1.3 and y = 
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1.1 values are: 

Emin — [ufsmin (re^ r l) n ^2]^ (34a) 

EmaT = [Ufs max (re/ri)” 1 ^] 2 (34b) 

The normalized mass of the atmosphere blown-off, me/m a versus both Emin and arc 
shown in Fig. 13. 

Table 1 shows that the energy for atmospheric loss above a tangent plane is a small 
fraction (~1(H to 1(H) of the energy required to drive-off the entire atmosphere. For the 
Earth, this total loss energy is -10^ ergs and would be achieved via an impact of a lunar- 
sized object at 20 km/sec. In the case of Venus, the impact of a smaller -800 km radius 
object, at -20 km/sec, will drive off the atmosphere. For Mars, the impact of a 160 km 
radius object at 20 km/sec will drive off the atmosphere. It may be, as suggested by 
Cameron [1983], that the terrestrial planets all suffered several giant impacts and then- 
present atmospheres may reflect accretion and outgassing since the last great impact event. 
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FIGURE CAPTIONS 

Figure 1. Formation of asteroid-size intermediate bodies form the dust component of the 
solar nebula a conditional variant assuming (a) instantaneous formation of the solar 
nebula and (b) breakup of the dust disc caused by the onset of gravitational 
instability. There is gradual accretion of intermediate bodies into planets. The 
accretion of gas by giant planets is not shown. The initially flat system of 
intermediate bodies thickens due to their mutual gravitational perturbations, (after 
Levin, [1972]). 

Figure 2. Enhanced image of the star (i-Pictoris demonstrates what appears to be the 

beginnings of another solar system. The disk of material surrounding {S-Pictoris 
extends 60 x 10 9 km from the star, which is located behind a circular occulting 
mask in the center of the image. The disc material is probably composed of gases 
and grains of ices, carbonaceous chondrite-like organic substances, and silicates. 
These are the materials from which the comets, asteroids, and planets of our own 
solar system are thought to have formed (after Smith & Terrile [1984]). 

Figure 3. Atomic abundance of the elements in the solar photosphere vs. the abundance in 
chondridc meteorites. Plot is normalized with respect to 10 6 atoms of Si (after 
Allfcgre, [1982]). 

Figure 4. Abundances of noble gases in planetary atmospheres and chondritic meteorites 
given as cubic centimeters per gram of rock (after E. Anders, personal 
communication in Owen et al [1992]). 
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Figure 5. Cartoon indicating the difference between the thermal regime of accretion via 
"thermal blanketing" and previous studies (figure from Abe & Matsui [1985]). 

Figure 6. (a) The evolution of surface temperature during accretion of a model Earth from 
planetesimals containing 0.1% H 2 O. The radius, R, is normalized by the final 
value, Rq. The dashed curve gives the calculated surface temperature without an 
impact-generated atmosphere (accretion period is 5 x 10 7 yr). The model surface 
temperature is affected by an atmosphere (which begins to greatly increase its mass 
once the impact velocity exceeds a critical value). The rapid rise in the surface 
temperature of the "standard model" which occurs after the Earth grows to -0.3 Ro, 
is due to an increase in the total mass of the atmosphere because of the initiation of 
the complete dehydration reaction of the surface layer. Once the surface 
temperature reaches the melting temperature, it remains nearly constant (figure after 
Matsui & Abe [1986a]). (b) The total mass of the impact-generated H 2 O 
atmosphere is plotted against the normalized radius for the standard planetesimal 
models. Note that the total atmospheric mass, M s , remains nearly constant after the 
Earth grows to 0.4Ro and is very close to the present mass of the Earth’s oceans 
(1.4 x lO^kg) (figure after Matsui & Abe [1986a]). 

Fig. 7. Shock front at successive instants of time for a strong explosion at high altitude. 
Sections shown are formed by passing a vertical plane through the origin of the 
explosion. The density of the atmosphere changes by a factor of e over the 
atmospheric scale height, A. Note that x = (p c AfyE) 1 ^ 7 where p c is the density at 
the altitude of the explosion and E is the energy (after Zel'dovich & Raizer [1966]). 

Fig. 8. Cartoons of atmospheric impact indicating how projectile momentum can be 
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transferred to a planetary atmosphere. A) The projectile enters the atmosphere, 
heating, compressing, and accelerating the atmospheric gases ahead of it B) Solid 
ejecta from a growing crater pass through the atmosphere, transferring some or all 
of their momentum to the atmosphere by drag. Only a small quantity of atmosphere 
is ejected via the mechanisms of A and B. C) The impact-generated vapor plume 
expands upward and outward (after Vickery & Melosh [1990]). 

Fig. 9. The impact-generated gas is assumed to interact only with the atmosphere lying 

above the plane tangent to the earth at the center of impact (after Vickery & Melosh 
[1990]). 

Fig. 10. The normalized maximum atmospheric mass that can be expelled by an impacting 
veneer of normalized mass for three values of q, the power law exponent describing 
the (differential) mass spectrum of impactors. Mars (dots) and Venus (dashes) are 
quantitatively similar (both shown for q = 1.7). Rectangles indicate veneer masses 
(width corresponds to plausible q values) needed to remove present atmospheres of 
Titan, Mars, and Venus (after Zahnle et al [1988]). 

Fig. 11. (a) Sketch of lunar-sized planetesimal impacting the Earth. The protoatmosphere 
is blown away by the shock wave-induced motion of the solid or molten planet (b) 
Impact of asteroidal-sized impactor and resultant partial eroded atmosphere. 

Fig. 12. Relationship of outward free-surface velocity (uf S ) to outward atmospheric 

velocity (v e ) for a polytropic atmosphere with various values of y. (a) Earth, (b) 
Venus, (c) Mars. 

Fig. 13. Calculated fraction of atmosphere blown-off versus impactor energy for Earth, 
Venus, and Mars. Lower and higher energy curves for each planet corresponds to 
assumed poly tropic exponent of ideal gas of y= 1.1 and 1.3, respectively. 
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Table 1. Characteristics of Impact Erosion Models 
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Table 2. Terrestrial Planet and Atmosphere Properties 
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